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Using the Wilson formulation of lattice gauge theories, a gauge invariant grid discretization of 
a one-particle Hamiltonian in the presence of an external electromagnetic field is proposed. This 
Hamiltonian is compared both with that obtained by a straightforward discretization of the contin- 
uous Hamiltonian by means of balanced difference methods, and with a tight-binding Hamiltonian. 
The proposed Hamiltonian and the balanced difference one are used to compute the energy spectrum 
of a charged particle in a two-dimensional parabolic potential in the presence of a perpendicular, 
constant magnetic field. With this example we point out how a "naive" discretization gives rise to 
Q\ ' an explicit breaking of the gauge invariance and to large errors in the computed eigenvalues and 

ON ■ corresponding probability densities; in particular, the error on the eigenf unctions may lead to very 

poor estimates of the mean values of some relevant physical quantities on the corresponding states. 
On the contrary, the proposed discretized Hamiltonian allows a reliable computation of both the 
r-{ • energy spectrum and the probability densities. 
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PACS number (s): 73.20Dx, 72.15Rn 



>. . I. INTRODUCTION 

OO 

In the last few years, there has been increasing activity concerning numerical simulations of low energy quantum 
\ systems and in this context particular attention has been paid to mesoscopic systems such as low-dimensional semi- 
conductor nanostructures. Of special interest are "quasi zero" dimensional systems; in this case, the electrons are 
confined in a small region of space -as in an atom- and hence they exhibit a quantized energy spectrum. In particular, 
ON ■ urged by many experimental resultstl, theoretical investigations of the influence of a magnetic field on the spectrum 
' of such systems of confined electrons have been carried outo. n r-, 

The problem of computing the energy spectrum of such systems has been addressed in many waysaH. In most 
of the approaches used the resulting Schrodinger equation is solved numerically, introducing a discretization grid 
and using finite difference methods (see e.gJa). In order to study the effect due to the presence of a magnetic field, 
the common procedure is to introduce a vector potential ("gauge field") and eventually construct an appropriate 
discretization of the corresponding Hamiltonian. On the other hand, it is well known from lattice quantum field 
theory that the presence of a gauge field may lead to the breaking of the gauge invariance if the Hamiltonian is 
naively discretized on a grid. In particular the eigenvalues of the Hamiltonian and the probability density associated 
with a given cigenfunction (i.e. its square modulus) could depend on the gauge choice. 

The aim of this paper is to present a method of constructing a grid-discretized one-particle Hamiltonian in the 
presence of an external electromagnetic field, that leads to a manifest gauge invariance of the mean values, of the 
observables on its eigenstatps. In order to do so we apply the Wilson formulation of lattice gauge theoriesO (for a 
review on this topic see alsc£3). 

The paper is organized as follows: in Sec. H, we construct explicitly a lattice discretization of the one particle 
Hamiltonian in the presence of an external electromagnetic field using the above mentioned formulation and we 
compare it with that obtained by introducing directly the vector potential and with a tight-binding one. In particular, 
we write the proposed Hamilto nia n in the tight-binding formalism and we show that it leads to a simpler expression for 



the hopping potential. In Sec. Ill we use the "naive" discretized Hamiltonian and the gauge-invariant one to compute 
the energy spectrum of a charged particle in a two-dimensional parabolic potential in the presence of a perpendicular, 
constant magnetic field , showing explicitly that the Wilson formulation preserves the gauge invariance. Finally, 
Sec. contains some concluding remarks. 



II. GAUGE INVARIANT LATTICE HAMILTONIAN 



In this section we apply the Wilson formulation to find a Hamiltonian -Hlatt defined on a lattice of discretization 
points that in the continuum limit (i.e. lattice spacing going to zero) tends to the following generic Hamiltonian 
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operator for a particle with charge q in the presence of a magnetic field B(f ) = V A A(r) and of an electrostatic 
potential V(f ): 



-nCONT — - — 
Am 



2 +V(r), (1) 



-V 2 + ig{A(r),v} + g 2 A(f) 



where g = q/h, and the symbol {A{f ), V} stands for the sum of the anticommutators for the components of A and 
V. Throughout this paper, for simplicity we neglect the spin. 

The essential requirement we impose is that the lattice Hamiltonian TJlatt maintains all the properties of the 
continuum one for a gauge transformation of the vector potential. 

It is well known from elementary quantum mechanics that when the vector potential transforms as 

A->A + VA(r) (2) 

a generic wave function ^(r ), corresponding to a physical state of the system, transforms as 

*(f ) -> G(f )*(f ) , with G{r ) = exp{igA(f )} , (3) 

leaving unchanged the mean values of the observables on 'J. In the following, we shall refer to the above mentioned 
property as gauge invariance. More formally, under the gauge transformation (^|) the Hamiltonian Hqont behaves as 

ffcoNT -» G(f ) ffcoNT Gt (r ) , (4) 

and, as a wave function 'J transforms according to (||), the energy functional 

£ c [*]= [ drV\r)H C oNTW) (5) 

is left unchanged. 

We now introduce a uniform discretization grid with lattice constant a, such that a generic point on it is identified 
with a vector x — (la,ma,na); we shall also indicate a generic direction in the discretized space with the Greek 
letter /i . Given a vector potential with components A^{x), followingcl we define on the lattice the operators U^{x) as 
follows: 

U^(x) = exp{iga A^(x)} . (6) 
Let us consider the discrete gauge transformation corresponding to (||), i.e. 

A^ -> A^ + - [A(x + m) - A(a)] , (7) 
a 

where x + fi denotes the next neighbor of x in the /z direction. It is straightforward to show that under gauge 
transformations (^) U ^ behaves as 

U^(x) -> G f {x) Un(x) G(x + jtt) , with G(x) = exp{igA(x)} . (8) 

In the following, we represent the discretized Hamiltonian operator Hlatt with its matrix elements between the 
functions x( x ~ x i) = ^-/*J~Q?&x,Xi-> where 5 x ^ Xi is equal to zero anywhere on the grid except at the lattice point xi 
where it equals 1, and l/^fa? is a normalization factor. Hence we define ^?latt {x, y) as: 

H hATT (x, y) = a 3 ^ x f ~ x)#lattx(S - y) = ^ ^.s^LATT&s,^ (9) 

X X 

With the definition above we have: 

H hATT ^{x) =^2H hATT (x,y)y{y). (10) 

y 

We now want to construct a lattice Hamiltonian which transforms as 

#latt(z, y) - G(x)H LATT (x, y)G\y) . (11) 
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As in the continuum case, this implies that a wave function \l/(r ) transforms according to the discrete analog of (||): 

*(x) G(x)V(x) . (12) 

We find that the following Hamiltonian 



h 2 1 

%H(l, V) = 7T~ ^ ~ [ 2 ^-V ~ U ti(v) S X-ll, V ~ U^X) S x+Ilt y] + V(x)5 x> y 



2m e — ' a 



(13) 



has the correct continuum limit (Jl|) and it transforms as required by Eq. ([11]). Moreover the Hamiltonian (|13|) gives 
rise, as in the continuum case, to a gauge invariant energy functional 



^latt[^] =a 3 ^ft^ffwHti,!/)*!!/)- 



(14) 



In order to show the advantages of the above formulation in respect to a straightforward discretization of (jl]) we now 
consider the discrete Hamiltonian obtained from (Q) applying balanced difference methods, which reads: 



H, 



BDH 



fj- 



(1 + iga Sx-^y - (1 - iga A^(x)) 8 x +n, y \ + V{x)6 XtV 



(15) 



One of the inconveniences of the Hamiltonian (|15|) is that it is not Hermitian and becomes such only in the continuum 
limit. 

As a discretized Hamiltonian is generally used to compute the lowest bound states of a system, we discuss briefly 
the error introduced by the discretization in this kind of computation. In particular we focus our attention on the 
error made when a discretized Hamiltonian acts on the eigenfunctions corresponding to the bound states we wish to 
compute. For this purpose we consider the quantity: 



A n (x) = 



*l(x)H C (mT(x)®n(x) - *n(«)-ffLATT(a?,y)*n(») 



(16) 



where ^ n is the exact wave function corresponding to the n-th eigenvalue of -ffcONT and all quantities are evaluated 
at a generic grid point x. It is worth noticing that the error A„ has a component due to the discretization of the 
Laplacian operator on the grid and an other due to the discretization of the gauge field. The systematic error due to 
the discretization of the spatial derivatives is 



A°(z) 



94 



2m* 



(17) 



This error is the same for both Hamiltonians ( |i"3| ) and (|l5|). As far as it concerns the error due to the discretization 
of the gauge field, the following considerations apply: When in Eq. ( pf ) for i?LATT we use i?wH (given by Eq. [l3| ), 
both ^(af)ifoONT(af) 1 S r n (a;) and ^2 y ^ (x)Hwh(x, y)^5>(y) are gauge invariant; hence the contribution to A„ due to 
the vector potential is gauge invariant too and it depends only upon the coupling between the particle and the field 
(i.e. g), upon the magnetic field B and upon the lattice spacing a. It may be shown that the component of the error 
A n due to the vector potential is small for any lattice point x if: 



\ga 2 B{x)\ « 1 



\/x . 



(18) 



Therefore, when condition ( |18| ) is satisfied and the systematic error ( |17| ) is small for the states we wish to compute, 
the Hamiltonian ([l3]) is a good approximation of ([!]). We notice that Eq. (|l^) is consistent with the fact that the 
stronger the magnetic field the more localized arc the eigenfunctions and thus a smaller lattice spacing a is required 
to get a good description of them. Instead, if in Eq. ( |l8| ) we use -ffsDH, then the quantity ^2 y (x)Hbdh(x, y)^{y) 

is no more gauge invariant and as a consequence the error A„ depends on the choice of the vector potential A. In 
particular for some "unhappy" gauge choices the error given by ( [l6| ) may be so large that (15) becomes useless even 
for a rough estimate of the ground state energy. 
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As tight-binding Hamiltonians are often used to discretize the Schrodinger equation on a lattice, it is interesting 
to see hpw the proposed Hamiltonian (|l3| ) compares to them. Without a magnetic field, a tight binding Hamiltonian 
(see e.gilJ) is equivalent to that obtained with balanced difference methods, although-written with another formalism. 
In the presence of a gauge field, a phase shift is introduced in the hopping potentiafij. A tight-binding Hamiltonian, 
for discretizing a Schrodinger problem in the presence of a magnetic field, reads: 



Htb — 2^ 



V(x) 



\x)(x\ 



EE 

X y^Lx 



2ma 2 J 



*){y\, 



(19) 



where Nd is the number of dimensions, the summation over y extends only to the next-neighbor of the lattice point 
x, and the phase shift is given by: 



A ■ dfi. 



(20) 



with the integration being performed on the straight line connecting the two points x and y. The first thing to 
notice is that Hamiltonian (M) has the correct continuum limit and is gauge invariant, as can be easily proved. The 
expression (^o|) for the phase shift was introduced for the tight-binding computation of electronic states in crystal 
lattices. For a crystal, the lattice constant is fixed and depends on the particular structure being investigated; therefore 
the gauge field A can vary considerably over a lattice constant and thus the use of an integral in the phase shift (|2p| ) 
is needed. This consideration does not apply to the case of a lattice of discretization points, as the lattice constant 
is chosen sufficiently small in order to get a good approximation. We also wish to point out that the proposed 
Hamiltonian (J13|) can be written in the tight-binding formalism: 



TBW 



E 



2ma 2 



V{x) 



k)(z| + E 



X ,/J 



2ma 2 



U^{x - n)\x){x - 



U fl \x)\x)(x + fi\ 



(21) 



Equation (|l|) suggests a simpler expression for the hopping potential in the presence of a magnetic field. In fact the 
magnetic field can be taken into account by introducing the lattice operator U^. Moreover, the Hamiltonian (|2l]) does 
not require the knowledge of the gauge field over the entire space but just on the lattice points; this is very useful 
if the vector potential is computed by some discretization technique. Finally, it is worth noticing that, as far as the 
discretization of a Schrodinger problem is concerned, Hamiltonian ( pl| ) is practically coincident with ([l9]), although 
Eq. ( |l9| ) is more complicated. 



III. NUMERICAL EXAMPLE 



In this section, we use Hamiltonians (13) and ( |15[ ) to compute the lowest part of the energy spectrum and the 
corresponding eigenfunctions for an exactly solvable problem: an electron confined in a two-dimensional parabolic 
potential in GaAs in the presence of a constant perpendicular magnetic field of strength B. The fact that this problem 
is analytically solvable allows us to test the accuracy of the discretized Hamiltonians and their behavior under gauge 
transformations. 

The Hamiltonian of the system is: 



H, 



CONT 



2m* 



{1(f), v} 
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A(f) 



1 



*n 2 i 



rc\ 



(22) 



where g = —e/h with e modulus of the electron charge, hfl is the confining energy, fc is the coordinate of the center 
of the parabolic potential, and m* is the effective mass for an electron in GaAs. 

To approach analytically the problem of computing the eigenvalues and the eigenvectors of the Hamiltonian (|2S 
it is convenient to use spherical coordinates on the plane where the electron is confined with the origin at the center 
of the parabolic potential. 

With the gauge choice A = — (l/2)r A B, the eigenvalues and eigenvectors of (E2n are given bylla: 



E n ,M = 2^ e ff^ n 



\M\ 



1) + -Uu) c M 



(23) 



where oj c = eB/m*, Q e ff = 4f2 + lj c , A = h/(m*fl e ff), n is a non negative integer quantum number, M is the 

angular momentum quantum number and L„' M ' is the n-th Laguerre polynomial. 

In order to compute numerically the lowest eigenvalues and eigenvectors of ( p2|) with the discretized Hamiltonians 
( |l3| ) and (|l5|), we introduce a finite square grid centered around the point fc, and we impose Dirichlet boundary 
conditions. This is equivalent to approximating the parabolic potential V with a potential V which is equal to V 
inside the grid region and goes to infinity outside it. We notice that the larger the region covered by the grid the more 
eigenvalues of (g2|) can be effectively approximated. In the following part of this section, dealing with the discretized 
computation, we use an orthogonal coordinate system with the grid placed in the first quadrant of the xy plane and 
the z-axis perpendicular to it. In all the numerical simulations we have taken a confining energy h£l of 4 meV and 
we have used for m* the value 0.067mo with mo being the electron mass. We have used a grid of 99 x 99 points and 
a lattice constant a of 2.4 nm, thus covering a region of 240 nm x240 nm. 

As our purpose is to investigate the behavior of the discretized Hamiltonians (|l^) and ( p"5| ) under gauge transfor- 
mations, we compute the eigenvalues and the eigenvectors for the following three gauge choices: 



(25) 



auf 


;e 1: 


A x 


= - V -B 
2 






Ay 


X 








auf 


;e 2: 


A.,. 


y-yc B 

2 






Ay 


= X -'°B 
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auj 


;e 3: 


A x 


= -(x-y)B 






A x 


= yB. 



(26) 
(27) 

These gauges are transverse; in particular, gauge 1 corresponds to A = — (l/2)r A B, while gauge 2 corresponds to 
A = — (1/2) (f—fc) A B. Gauge 2 is similar to gauge 1, but symmetric with respect to the center of parabolic potential 
and, as we shall see below, it yields more accurate results with respect to the other two gauges when inserted in the 
naively discretized Hamiltonian ([l5]). In the following, we shall refer to the Hamiltonians ( |l3| ) and ( p"5| ) as WH 
(Wilson Hamiltonian) and BDH (Balanced Difference Hamiltonian) respectively. 

In Fig. 1 we show the lowest fifteen eigenenergies for a magnetic field strength of 10 T, computed with WH and 
BDH, together with the exact ones. We have numbered the eigenvalues with an integer i so that Ei < Ei + i and 



we plot them versus i. We have found that the results obtained using WH with the three gauges (25 2q , p7| ) do not 
differ; therefore for the quantities computed with WH we do not specify which gauge we have used. From the above 
mentioned figure, we notice that: (1) The eigenvalues computed with BDH using gauges 1 and 3 differ appreciably 
from the exact ones; in particular, this difference grows with increasing eigenvalue. (2) The eigenvalues computed 
with BDH using gauge 2 tend to differ from the the exact ones for the excited states. 

The explicit breaking of gauge invariance present in the Hamiltonian can be further investigated analyzing the 
probability density I'I'^mI 2 : in Fig. 2 we show the probability density for the ground state (with B = 10 T) using 
either the corresponding exact eigenfunction (from Eq. p4| ) or that obtained with BDH; for this computation we have 
not found significant differences between the exact probability density and cither that computed with WH (using the 
three different gauges) or that computed with BDH using gauge 2. From Fig 2 it is evident that in this case the 
breaking of gauge-invariance of BDH "disrupts" the symmetry of the probability density with respect to the center 
of the confining potential. We have found the same phenomenon with the gauge 3, while this does not occur using 
gauge 2 because this gauge preserves the above mentioned symmetry. 

In Fig. 3 we report the probability density for the 14-th excited state, i.e. that with n = 0, M = —14 (see Eq. p3[ ), 
and B — 10T using the corresponding exact eigenfunction and that computed with BDH using gauge 2; also in this 
case we have not found significant differences between the exact calculation and the one performed with WH, while 
it is evident that the result obtained with BDH using gauge 2 does not match the exact one. At this point it is worth 
stressing that even if in some cases, for a careful gauge choice, BDH gives a good estimate of the eigenenergies, the 
estimate of the corresponding wave functions may be inadequate (see Fig 1 and Fig. 3). This implies that the matrix 
elements of an operator corresponding to an observable physical quantity (e.g. dipole and quadrupole momenta) 
between the eigenf unctions computed with BDH can be affected by large errors. 

Finally, in Fig. 4 we show the ground state energy as a function of the magnetic field magnitude computed with 
the exact expression and with the three different gauges using both WH and BDH. The results obtained using WH 
with the three gauges (^5 26 127]) do not differ; therefore for the quantities computed with WH the gauge used is not 
specified. This analysis shows how the effect of the breaking of the gauge invariance in the spectrum of BDH grows 
with increasing magnetic field magnitude. 
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IV. CONCLUSIONS 



In this paper we have addressed the problem of developing a method to find a gauge invariant grid discretization 
of a one-particle Hamiltonian in the presence of an external electromagnetic field; for this purpose, we have used the 
Wilson formulation of lattice gauge theories. This Hamiltonian, written in the tight-binding formalism, leads to a 
simpler expression for the hopping potential, that does not require the knowledge of the gauge field over the entire 
space but only on the lattice of discretization points. The method developed has been checked in an exact solvable 
case, i.e. a Hamiltonian describing an electron confined in a two-dimensional parabolic potential in the presence of 
a constant perpendicular magnetic field. The results obtained clearly show how a "naive" discretization yields to 
an explicit breaking of the gauge invariance and may lead to large errors in the estimates of the eigenvalues and 
corresponding probability densities. With a careful gauge choice a naively discretized Hamiltonian can be used to give 
a rough estimate of the energy spectrum, but the corresponding eigenfunctions are still affected by large errors, thus 
being unsuitable for the computation of the matrix elements of operators corresponding to relevant physical quantities. 
On the other hand, the proposed discretized Hamiltonian ( |l3| ) has the correct behavior under gauge transformations 
and allows an accurate computation of both the energy spectrum and the probability densities. 
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FIG. 1. Lowest fifteen eigenvalues for B = 10T computed with: exact Hamiltonian (points); WH (diamonds); BDH with 
gauge 1 (circles); BDH with gauge 2 (cross); BDH with gauge 3 (squares). 



FIG. 2. Ground state probability density for B — 10T: a) exact; b) computed with BDH (gauge 1). 



FIG. 3. Probability density of the 14-th excited state (n = 0, M = —14): a) exact; b) computed with BDH (gauge 2). 



FIG. 4. Ground state energy plotted against magnetic field magnitude computed with: exact Hamiltonian (continuous line); 
WH (diamonds); BDH with gauge 1 (circles); BDH with gauge 2 (cross); BDH with gauge 3 (squares). 
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